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ABSTRACT 

We report the measurement of the primordial D/H abundance ratio towards 
QSO HS 0105+1619. The column density of the neutral hydrogen in the 
z ':=L 2.536 Lyman limit system is high, log A"/// = 19.422 ± 0.009 cm~^, 
allowing for the deuterium to be seen in 5 Lyman series transitions. The 
measured value of the D/H ratio towards QSO HS 0105+1619 is found to be 
D/H= 2.54 ± 0.23 x 10^^. The metallicity of the system showing D/H is found 
to be ~ 0.01 solar, indicating that the measured D/H is the primordial D/H 
within the measurement errors. The gas which shows D/H is neutral, unlike 
previous D/H systems which were more highly ionized. Thus, the determination 
of the D/H ratio becomes more secure since we are measuring it in different 
astrophysical environments, but the error is larger because we now see more 
dispersion between measurements. Combined with prior measurements of 
D/H, the best D/H ratio is now D/H= 3.0 ± 0.4 x 10-^ which is 10% lower 
than the previous value. The new values for the baryon-to-photon ratio, and 
baryonic matter density derived from D/H are 77 = 5.6 ± 0.5 x 10^^° and 
^]6/^2= 0.0205 ± 0.0018 respectively. 
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Subject headings: quasars: absorption lines - quasars: individual 
(HS 0105+1619) - cosmology: observations 

1. Introduction 

The standard theory of big bang nucleosynthesis (SBBN) predicts the 
abundances of the hght nuclei H, D, ^He, '^He, and ^Li as a function of the 
cosmological baryon to photon ratio, rj — nb/riy (Kolb & Turner 1990; Walker 
et al. 1991; Schramm & Turner 1998; NoUett & Buries 2000; Olive, Steigman 
& Walker 2000). A measurement of the ratio of any two primordial abimdances 
gives 77, and hence the baryon density, while a second ratio tests the theory. 
However, it is extremely difficult to measure primordial abundances, because in 
most places gas ejected from stars has changed the abundances. 

Adams (1976) suggested that it might be possible to measure the primordial 
D/H ratio in absorption hue systems towards QSOs. Although gas which has 
been inside a star will have lost all of its deuterium, in QSO absorption line 
systems having typical metal abundances of 0.001 to 0.01 of the solar value, 
about 0.1 - 1% of the deuterium will have been lost. 

The advent of the HIRES spectrograph (Vogt 1994) on the W.M. Keck-I 
telescope gave the high signal-to-noise ratio (SNR) and spectral resolution 
needed to reveal deuterium (Tytler ct al. 2000) in high rcdshift absorption 
systems. We have previously measured D/H in two QSOs (Tytler, Fan & Buries 
1996; Tytler & Buries 1997; Buries & Tytler 1998a; Buries & Tytler 1998b), and 
placed a strong upper limit on D/H in a third (Kirkman et al. 1999). 

Other QSOs give less useful constraints on D/H, because their absorption 
systems are more complex, or existing spectra are inadequate. The Lyman limit 
system (LLS) at Zabs = 0.701 towards QSO PG 1718+4807 might allow ten times 
larger D/H, or it may give no useful constraints (Webb et al. 1997; Levshakov, 
Kegel, & Takahara 1998; Tytler et al. 1999). Molaro et al. (1999) claimed 
another QSO absorption system showed low D/H, but they and Levshakov et 
al. (2000) note that since only the Lya line has been observed, the hydrogen 
velocity structure and the H I column density are poorly known, and the 
deuterium feature can be fit using hydrogen alone. 

In this paper, we present a fourth QSO, HS 0105+1619, which gives strong 
constraints on the primordial D/H ratio. 
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2. Observations and data reduction 

We report the detection of deuterium in tlie QSO HS 0105+1619 (emission 
line redshift 2.64, V=16.9, B1950 RA Ih 5m 26.97s, DEC +16d 19m 50.1s; 
J2000 RA Ih 8m 6.4s, DEC +16d 35m 50.0s), which was discovered by Hagen, 
Engels & Reimers (1999), who very kindly gave us a finding chart and a low 
resolution spectrum prior to publication. 

We present high quality, high SNR spectra of HS 0105+1619 in both low 
and high resolution. The low resolution spectra were obtained using the Kast 
double spectrograph on the Shane 3 meter telescope at Lick observatory. The 
high resolution spectra were obtained using the HIRES spectrograph on the 
Keck-I telescope. The observations are summarized in Table |l]. All of the high 
resolution observations were taken using the HIRES image rotator to align 
the direction of atmospheric dispersion along the slit, and were taken with 
the C5 decker, which provides an entrance aperture to the spectrograph with 
dimensions 1.15" x 7.5" . The spectra were sampled in 2.1 km s~^ pixels with 
the Tektronix 2048x2048 CCD. 

The HIRES spectra were fiat-fielded, optimally extracted, and wavelength 
calibrated using Tom Barlow's set of echelle extraction MAKEE programs. The 
spectra were then co-added to produce a single spectrum. This spectrum was 
then fiux calibrated using the low resolution Kast spectrum. Because the SNR 
decreases in both spectra at lower wavelengths, we use the fiux calibrated data 
for wavelengths greater than 3800 A and the co-added spectrum without fiux 
calibration below 3800 A. The details regarding the co-adding and fiuxing are 
given in Suzuki & Tytler (2000). The wavelength scale has an external zero point 
error of approximately ±10 km s~^, which is the shift between spectra taken at 
different times, which we corrected. The internal error in the wavelength scale 
is a least 0.09 km s~^, from the arc line fits, and may be approximately 1-2 
km s~^, which is the size of the error which we have noted in our analysis of 
other similar spectra (Levshakov, Tytler & Buries 2000). The final spectrum 
has 2.1 km s~^ wide pixels, and a SNR of approximately 80 and 10 at the Lya 
and Lyman limit of the D/H system, respectively. In Figure 0, we show the fiux 
calibrated Kast and HIRES spectra. 

3. General Properties of the z ~ 2.536 Lyman Limit System Towards 

HS 0105+1619 

As can be seen in Figure |I|, the low resolution Kast spectrum of HS 0105+1619 
shows a steep Lyman break at a wavelength of approximately 3230 A, which is 
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caused by a Lyman limit system at a redshift of z ^ 2.54. At this redshift we see 
a strong Lya absorption line, near 4300 A. Our analysis of the high resolution 
HIRES spectrum indicates that there is a Lyman limit system at a redshift of 
z ~ 2.536, which shows deuterium and numerous metal absorption lines. We 
discuss the general features of the absorption below, and the best fits to the 
data in the following section. 

3.1. Hydrogen Absorption 

We observe hydrogen in all Lyman series transitions through Lyman- 17, 
where the spectrum abruptly ends at 3233 A due to line blending near the 
Lyman limit. The observed Lyman transitions are shown in Figure ^ 

The column density of the hydrogen is high, as Figure |^ indicates. All of the 
Lyman series transitions are saturated through to the Lyman limit, indicating 
that the column density is at least log N^i > 17.8 cm~^. The Lya line has zero 
flux across about 200 km s~^ indicating that the column density of the system 
is approaching the levels found in damped Lya systems. 

Inspection of Figure || also indicates that the absorption system is very 
simple. The Ly-5, Ly-6, Ly-7, Ly-9, Ly-10, Ly-14, and Ly-15 transitions all 
appear symmetric and un-blended, allowing us to describe the absorber by a 
single component. 

3.2. Deuterium Absorption 

Since the column density of the hydrogen responsible for the Lyman limit 
appears high, we expect the associated deuterium absorption to be strong. For 
the first time, we observe absorption at the predicted position of deuterium, 
a velocity of f = —81.64 km s~^ in the frame of the H I, in 5 Lyman series 
transitions: Ly/5, Ly7, Ly-5, Ly-6, and Ly-7. Deuterium Lya was not observed 
since it is subsumed by the absorption of the hydrogen Lya, and deuterium Ly-4 
was not observed due to intervening Lya forest absorption. 

The observed absorption is narrow, and appears symmetric and free of strong 
contamination, suggesting that like the hydrogen, the absorption is simple and 
can be modeled with a single component. 

In all transitions where absorption is observed, the velocity of the absorption 
appears centered about the same value, v ~ —82 km s~^ , strongly indicating 
that the features are all Lyman transitions of the same absorber, and that the 
absorber is consistent with being deuterium. In a later section, we give a more 
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rigorous discussion of why we believe the absorption is indeed deuterium, and 
not hydrogen or other ions. 

3.3. Metal Line Absorption 

The z ^ 2.536 Lyman limit system shows a variety of metal ions, as seen 
in Figure ^ Unlike Lyman Limit systems with column densities in the range 
of \o'gNHi= 16.5-18 cm^^ which show absorption predominantly in the higher 
ionization states, this system shows metal line absorption in neutral, low, and 
high ionization states. Since the column density of the hydrogen appears high, 
we expect the neutral and low ionization metal lines to trace the H I, and with 
it the deuterium, in analogy with damped Lya systems. 

All of the low ionization metal lines are extremely narrow, and appear to be 
well described by a single component. Moreover, the low ions are all centered at 
~ km s~^, implying that they arise in the same gas as the H I and D I. 

4. Best Parameters for the z 2.536 Lyman Limit System 

We now give the parameters which describe the absorption in the z ~ 2.536 
Lyman Limit system towards HS 0105+1619. For all measurements, a continuum 
was fit to the region under consideration to produce a unit normalized spectrum. 
The observed absorption features were then fit using the VPFIT Voigt profile 
line-fitting routine (Webb, 1987; kindly provided by Carswell) and re-verified 
using in-house routines. For each absorber, we obtain an estimate of the column 
density. A, the redshift, 2, and the velocity width, 6, along with their respective 
lex errors. The results of this analysis are found in Table 

4.1. The Hydrogen 

The parameters describing the hydrogen absorption responsible for the 
Lyman limit are obtained from various complementary parts of the spectrum, 
and are given by log Ah/= 19.422 ± 0.009 cm'^, h = 13.99 ± 0.20 km s-\ and 
z = 2.535998 ±0.000007. 

The deuterium and low ionization metal absorption lines give a strong 
indication that the H I can be modeled by a single component, whose redshift 
should be consistent with that of the other neutral ions observed, namely O I 
and N I. We find that a single component fit allows for an excellent description 
of the observed absorption, and we now consider the individual parameters. 
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4.1.1. Redshift 

The redshift of the hydrogen was determined by simultaneously fitting the 
Lyman series absorption at Lya, Ly/?, Ly7, Ly-5, Ly-6, Ly-7, Ly-9, Ly-10, 
Ly-13, Ly-14, and Ly-15. 

4.1.2. Velocity Width 

The high column density of the absorber does not give us information on 
the velocity width at Lya, so we must turn to the higher order Lyman series 
lines. To determine b, we fit the absorption in the Lyman series in transitions 
which appear to be least contaminated by interloping Lja forest, specifically 
Ly-9, Ly-10, Ly-14, and Ly-15, simultaneously. The inclusion of all observed 
Lyman series transitions in the fitting procedure does not changed the measured 
b value. 

4.1.3. Column Density 

Three regions of the spectrum allow us to determine the column density: the 
Lyman limit, the core of the Lya, and the damping wings of the Lya. These 
regions are shown in Figure § along with the single component best fit to the 
H I. 

The Lyman limit allows for a lower limit to the column density of the 
hydrogen of \ogNHi= 17.8 cm~^ since all observed Lyman series lines observed 
are saturated. 

The Lya line gives the most information regarding the hydrogen column 
density, because it is insensitive to the b value, even though it is sensitive to 
the continuum level. In Figure ^, we show the continuum used to fit the 250 
A region encompassing Lya and the corresponding approximate Icr levels, 
which amount to a ±2% continuum level change. To determine the effect of 
the continuum placement on the measured value on the column density, for all 
regions of the Lya line, we obtain a fit using the best continuum estimate, and 
then move the continuum to the itlo" levels and re-fit. In all fits to the Lya, we 
choose segments of the spectra which appear to be least contaminated by other 
absorption. 

In Figure ^, we show the fit to the core of logA'^^/= 19.419 ± 0.009 cm~^, 
where the error is the quadratic sum of the error from the continuum (0.007) 
and the error from the fit (0.006). This fit was made for two regions on either 
side of the line center: 4294.5-4295.5 and 4300-4302 A. There is additional 
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absorption in the core between 4295.5 and 4297 A which is readily fit by two 
H I Lya hnes. Their absorption is seen and fit in Ly/?, but they contribute 
no significant optical depth to the lower wavelength (4294.5-4295.5 A ) region 
which gives the column density of the H I which shows deuterium. 

We note that contamination of the core region is possible, and would 
lower the measured column density. However, we consider such contamination 
unlikely, since it would require at least two lines appearing in the right places, 
with a very restricted set of parameters, to produce enough absorption to fit 
both sides of the core region. 

The damping wings of the Lya absorber also give the column density, but 
the exact continuum placement is now the dominant source of error because 
the continuum uncertainty represents a larger fraction of the total absorption 
in these regions. Two damping wing regions were fit, one on either side of 
the line center. On the blue (lower wavelength) side, we fit the wavelengths 
4288.4-4289.3 A, which gave log Nhi = 19.406 ± 0.060 cm'^ where the errors 
are the range in column density allowed by the ±2% range in the continuum 
placement. For the red (higher wavelength) side regions 4306-4307 and 
4308-4309 A were fit giving log Nhi = 19.4752 ± 0.040 cm'^. 

For our best value for the H I column density, we take the weighted mean 
of the three measurements from the core and wings: log Nhi= 19.422 ± 0.009 
cm~^. All three estimates of the column density are consistent with this mean. 

4.2. The Deuterium 

The determination of the parameters describing the deuterium absorption is 
relatively straightforward, given the presence of 5 un-blended transitions, most 
of which are unsaturated. We begin by fitting the D I transitions, and then 
assess possible sources of error. 

The D I transitions are well fit by an absorber with log Ndi= 14.810 ± 0.029 
cm~^ and a velocity width of 6 = 9.93 ± 0.29 km s~^. This fit is shown in 
Figure |^. These values were determined by fitting the regions surrounding 
the deuterium in all 5 deuterium absorption regions simultaneously, with the 
redshift of the deuterium tied to the hydrogen, and the parameters for the 
hydrogen absorption fixed at logA'^;^/= 19.42 cm~^ and b = 13.99 km s~^. In 
all fits, the regions used to fit deuterium were constrained to be those within 
approximately —150 and +20 km s~^ of the H I line center, and are listed in 
Table 0. Where needed, additional Lyo; forest absorption was fit to model all 
absorption. The parameters for the additional lines are given in table ^. 
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The uncertainty in the continuum level increases the error on the D I column 
density by a third. Independent estimates of the continuum levels at the D I 
transitions had a Icr dispersion of approximately 10% , much larger than the 2% 
error near Lya, where the data has been flux calibrated and has significantly 
higher SNR. To gauge the effect on the error in the D I column density, we 
determined the parameters of the deuterium absorption independently for 
each line, both for the best estimate of the continuum and for a 10% higher 
continuum level. The difference between these two is the contribution to the 
error from the choice of continuum level. For each each D I transition, this error 
was added in quadrature to the error obtained with the best continuum. The 
weighted mean for the five D I transitions gives log A'^d/= 14.826 it 0.039 cm~^ 
and b = 9.85 ± 0.42km s~^. The results of the independent fits to the different 
deuterium transitions with the best estimate of the continuum level can be seen 
in Figure ^. 

The error on log Ndj is insensitive to a change in the redshift of 1 km s^^, as 
found by both letting the deuterium fit freely to its own redshift, and by tying 
the redshift to the H at 1 km s^^ away from the best fit value. Both results gave 
a log Ndi which deviated by amounts significantly lower than the best fit la 
errors. 

The effect of varying the main hydrogen by the la errors in either column 
density or velocity width produced no effective change on the log N^j since the 
main hydrogen component is well separated from the deuterium because the b 
of the hydrogen is very small. 

4.3. The Metals 

The results of fits to the many metal lines are given in Table 0, and a subset 
of the fits is shown in Figure ^. In all cases, the ions were best fit by a single 
component whose redshift, velocity width, and column density were all allowed 
to vary. Many of these ions show multiple transitions, and for any single ion, all 
transitions present in the spectra were fit simultaneously. 

The metals were all found to lie within 8 km s~^ of the redshift of the 
hydrogen which shows deuterium. More importantly, since the column density 
of the hydrogen is so high, we expect the redshifts of the low and neutral ions to 
agree with the hydrogen since we expect the gas to be predominantly neutral. 
The neutral and singly ionized ions agree with the hydrogen redshift to within 
2 km s~^, while the neutral ions alone, in the top section of Table ^ agree 
to within approximately 1 km s~^. The larger dispersion seen for the singly 
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ionized ions, in the second section of the table, indicates that some of the gas 
making these hnes is distinct from the neutral gas, but this dispersion could be 
insignificant, because the internal wavelength errors could be 1 -2 km s~^. 

The detection of O I is of particular importance in this system for two 
reasons. First, the ionization potential of O I is nearly identical to that of 
hydrogen, so it is ionized by the same photons that ionize the H I. Second, O 1 
participates in electron transfer with H I, such that in cases where the gas is not 
highly ionized, O I/O is nearly identical to H I/H, and the distribution of O I 
should match that of the H I and the D I. The O I absorption gives information 
about the temperature, bulk motion, ionization and abundances in the neutral 
gas which shows the D I absorption. Since the O I is accurately modeled by a 
single, narrow component in four separate transitions, we gain confidence that 
a single component fit to the deuterium and to the hydrogen is sufficient. The 
implications of the measurement of O I on the metallicity and ionization are 
discussed in a later section. 

Since O I and H I should arise in the same gas, we use these lines 
to obtain the temperature of the gas which shows D 1 and its turbulent 
velocity. We model the observed 6 as 6^ = fe^^^^ + The instrument 
line broadening, 6?^^^ was measured from arc lamp calibration spectra to 
be binst = 4.81 ± 0.14 km s~^. We model the intrinsic velocity width as 
^Int — ^temp + ^turb'j ^ Combination of thermal broadening and bulk motion. The 
thermal broadening, b^^,^^ = 2kT/m = 166.41 (T/10^K)/mass(amu), depends on 
the ion mass in atomic units, m, but the bturb is the same for all ions. All of the 
b values quoted in this paper and listed in Table § refer to intrinsic line widths, 
but the listed errors do not include the error in binst, because we do not know 
whether this error is correlated at different wavelengths. 

Fitting O I, N I and H I alone gives T = 1.15 ± 0.02 x 10^ K, and 
bturb= 2.56 ± 0.12 km s~^, which we show by the straight line in Figure |T0|. The 
errors quoted here are very much minimum values, because they do not include 
the error in the binst or the appropriateness of the model. 

Ions C II, Si II and Fe II, are all wider than predicted by this fit, presumably 
because a part of each line arises in gas with different velocity structure. The 
lines from the higher ionization ions C III, C IV, N II, Si III and Si IV have 
velocities which differ by to -7 km s~^ from the H I and low ionization ions. 
They are not relevant to the D/H because their ionization and velocities show 
that they arise in different gas, a common finding for absorption systems with 
high Nhi (Wolfe & Prochaska 2000a), but unlike absorption systems with much 
lower Nhi, including PKS 1937-1009 and Q1009+2956. 
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5. Is the Observed Absorption Deuterium? 

Now that we have determined the parameters of the absorption at the 
position of deuterium, we turn to the issue of confirming that the absorption is 
indeed deuterium, and not inter-loping hydrogen or metal line contamination. 

The primary concern with the absorption seen at the position of deuterium 
would be that it is caused not by deuterium, but instead by inter-loping 
hydrogen. Here we argue that this scenario is unlikely for the following reasons. 

First, hydrogen lines in the Lya forest with the appropriate column density, 
log Nhi — 14.8 cm~^ have 6 > 20 km s~^, and not 6 ~ 10 km s~^ (Kim et al. 
1997; Kirkman & Tytler, 1997). However, such low values of b might be found 
in components of the LLS, which are the most likely contaminants. Second, 
we are able to predict the width of the deuterium using the measured widths 
of the other neutral ions observed which are present in the same gas. Figure 
|T0| illustrates the concept. Since we have observed three other neutral ions 
(H I, O I, N I), we can use their widths to predict the width of deuterium. As 
Figure |l^ shows, the measured value of b{D) is consistent with its predicted 
value. Third, hydrogen lines with a log Nhi— 14.8 cm~^ often show associated 
metal line absorption, but as is seen in Figure ^ there is no such absorption at 
—82 km s~^. Finally, for the absorption to be hydrogen and not deuterium, its 
position would have to agree with that of deuterium to within 1 km s~^. Taken 
together, these arguments indicate that the observed absorption is deuterium, 
and not hydrogen, but we can not quantify this because we do not know the 
properties of components of Lyman limit systems. 

The scenario in which the absorption is metal line contamination is even 
less likely for a number of reasons. Any metal lines with a column density 
of log Nmetai — 14.8 cm~^ would show absorption in not only that ion, but 
many others along with strong associated hydrogen absorption, some of which 
would be easily observed in our spectrum, but were not. Also, for the observed 
absorption to be entirely derived from metal line contamination, such metal 
lines would have to appear in 5 different regions of the spectrum, all at positions 
within approximately 1 km s~^ of the predicted positions of deuterium, and 
with line strength scaling as the oscillator strengths expected for the deuterium 
Lyman series. A similar argument can be used to exclude the case whereby 
the absorption at the position of deuterium was hydrogen, but not in the 
corresponding Lyman series transition (e.g., the absorption at Ly-6 is an 
unrelated Lya line). 
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6. Elemental Abundances and Ionization State of the z ^ 2.536 

Lyman Limit System 

Here we discuss the abundances of the elements observed in the z ~ 2.536 
Lyman limit system which shows D/H. 

The general procedure for determining the elemental abundances is to use 
an ionization model to convert from the observed column densities of selected 
ions into elemental abundances. On the whole, the level of ionization in QSO 
absorption line systems tends to decrease as the column density of the gas 
increases, since the gas can shield itself more from ionizing Lyman continuum 
radiation. Typical Lyman limit systems with log Nhi — 17.5 cm~^ are highly 
ionized, whereas damped Lya systems, whose column densities are greater than 
log Nhi = 20 cm~^ are typically neutral. In the case of HS 0105+1619, we 
expect the gas to be predominantly neutral (H 1/H> 0.5), since the column 
density is approaching that of a damped Lya system. 



6.1. Ionization &i Metallicity 

We modeled the ionization with the CLOUDY v90.04 package developed 
by G. Ferland (Ferland, 1991). As usual, we assumed a plane-parallel 
geometry with constant density uh = 0.01 cm~^, and we approximated an 
isotropic background ionizing radiation by placing a point source at a very 
large distance. We used the Haardt-Madau ionizing spectrum (Haardt & 
Madau, 1996). In Figure 11 we show the correction which would be added 
onto an observed ionic abundance to give an elemental abundance: [X/H] = 
log{Nx/NH i) + loge(X/H I) — log(A/if)0, where Nx is the column density of 
the ion used to infer the abundance of element X. Since we hold the Nhi, the 
gas density and the spectrum constant, these corrections are specified in terms 
of the intensity of the ionizing flux at 1 Ryd, log J912. Because the corrections 
are nearly homologous with the ratio nn/Jgu, we can account for different uh 
values by simply shifting J912. 

For the D/H absorber, we can obtain the [0/H] without ionization correction, 
because we have four un-blended O I transitions which give the most accurate 
Noi in any absorption system to date. The ionization correction for O I is 
negligible, loge(0 I/H I) < 0.2 dex, even at extreme values of J912, and O I/H I 
is equal to 0/H, which gives [0/H] = —2.0. 

This abundance is supported by four other elements. Figure |lT] shows that 
the ionization corrections for Si II, C II, Al II, and Fe II are all negative. 
Therefore, we can place upper limits on the C, Si, Fe, and Al abundances 
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using the observed C II/H I, Si II/H I, Fe II/H I and Al II/H I ratios without 
corrections. We find [C/H] < -1.9, [Al/H] < -2.1, [Fe/H]< -1.9, and 
[Si/H]< —1.85, all within 0.15 dex of [0/H]. Note that these elements might 
have shown different abundances, because they are seen in lines which have 
different velocities and velocity widths from the H I, D I and O I. From N I, 
we find [N/H] ~ —3.1 without ionization correction, as suggested by the 
other elements. This implies that N is under-abundant by approximately 1.1 
dex, which is not unusual for damped Lya systems with similarly low metal 
abundances (Centurion et al. 1998; Lu, Sargent & Barlow 1998), but is much 
more N than usual in Galactic and extra-galactic H II regions (Henry, Edmunds 
& Koppen 2000). 

The similarity of these abundance values shows that each element is primarily 
in the listed ion, O I, Si II, Fe II, Al II and C II, and that H I/H > 0.8. This 
low ionization requires {Jgu/lO^"^^ erg cm~^ s~^ Hz~^ sr"-*^) {10~'^ cm.~^ /uh) < 1, 
and hence a weak radiation intensity and high gas density. When we divide the 
Nhi by the gas density, we find that size of the cloud along the line of sight 
would be Z < 1 kpc, assuming a homogeneous medium. If the absorbing region 
were spherical, its gas mass would be < 1.2 x lO^M©. 

The observed Si III, C III, C IV or Si IV absorption must come from 
separate, more highly ionized gas. 

7. Best Fit Values for HS 0105+1619 

For the gas giving D/H in HS 0105+1619, the best values for various 
parameters, and the la errors, are as follows: 

• log Nhi= 19.422 ± 0.009 cm-^ (2% error) 

• log Ndi= 14.826 + 0.039 cm'^ (9% error) 

• log D/H = -4.596 + 0.040 (10% error) 

• D/H = 2.54 + 0.23 x 10"^ 

• temperature: T = 1.15 + 0.02 x 10^ K 

• Gaussian turbulent velocity width: bturb— 2.56 + 0.12 km s~^ 

• log oxygen abundance, on the solar scale: [0/H] = —2.0 

• other limits: [C/H] < -1.9, [Al/H] < -2.1, [Fe/H] < -1.9, [Si/H] < -1.85, 
[N/H]~ -3.1 

• neutral fraction: H I/H> 0.8 
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• gas density: n > 0.01 cm~^ (for J912 > 10~^^ erg cm~^ s~^ Hz^^ sr^^) 

• extent of absorbing region along the line of sight: < 1 kpc 

• mass of gas: < 1.2 x 10^ Mq. 

For D/H, we added the errors on log Nhi and log N^i in quadrature, because 
we expect little correlation. 

8. Best values for D/H from all QSOs 

The weighted mean D/H from three QSOs, PKS 1937-1009, Q1009+2956, 
and HS 0105+1619 (the other QSO, Q0130-4021, gives a consistent upper limit) 
is logD/H = —4.523 ± 0.026 (6% error). The individual parameters from our 
measurements of four QSOs are summarized in Table |^. The dispersion in the 
three measurements of D/H is larger than expected. We obtain a. X2 — "^-^^ 
which would be exceeded in 3% of random samples with Gaussian errors. There 
are two interpretations. Either we have under-estimated the errors, or there is 
real dispersion in D/H. We favor the former interpretation, since there are a 
number of systematic errors which affect each measurement, and are difficult to 
assess. For example, we may have under-estimated the effect of contamination in 
PKS 1937-1009 or Q1009+2956. Hydrogen absorption near the D I lines would 
artificially increase the D/H estimated for these objects. We now discuss another 
possible source of scatter in D/H, astration, the destruction of deuterium in 
stars. 

8.1. Primordial Nature of the Gas & Astration 

At a metallicity of 0.01 solar, all chemical evolution models which are 
consistent with Galactic data predict that the destruction of deuterium is about 
1%, ten times smaller than our measurement error (Tosi et al. 1998; Jedamzik 
& Fuller 1997). 

We can get an approximate estimate of the amount of deuterium which 
will have been destroyed in stars by linearly scaling from the local interstellar 
medium (LISM). Since 50% of the primordial deuterium is lost from the LISM 
where the metal abundance is approximately 0.7 solar, we expect 0.7% to be lost 
at 0.01 solar. This scaling is only approximate for many reasons: the relation 
is not linear in the LISM today because the destruction is saturating, there 
has been infall of less processed gas, and the high mass stars which we expect 
polluted the absorption systems destroy less deuterium for a given amount of 
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metal production. 

We now discuss the relation between the D/H values and the metallicity 
which we have measured towards QSOs. The D/H values for the four QSOs 
listed in Table grange from < 6.7 xlO~^ to 2.54 xlO~^, while the [Si/H] ranges 
from —2.6 to —1.85. Figure |12| shows that these values appear correlated, in the 
sense that the D/H is lower when the metallicity is higher. This is the sense 
expected if the deuterium has been destroyed in those absorbers with higher 
metalhcity. The D/H for Q1009+2956, PKS 1937-1009 and HS 0105+1619 
appears to decline linearly with increasing [Si/H]. The D/H declines by a factor 
of 0.64 between Q1009+2956 and HS 0105+1619 while the [Si/H] increases from 
approximately —2.53 to —1.85. 

There are two arguments against the reality of the correlation. First, we do 
not expect such a large destruction of deuterium when the metal abundance 
remains this low. A decrease in the D/H by a factor of 0.64 implies that 0.64 
of the atoms in the absorbing gas have not been in a star, and that 0.36 have 
been in one or more stars, which is at least an order of magnitude more than 
expected. Second, the indicated rate of decrease of the D/H is much faster than 
expected. We expect D/H to decrease linearly with the linear, and not log, 
metal abundance. 

Models of the evolution of the chemical abundances in our Galaxy predict 
that D/H will remain near the primordial value - a horizontal straight line on 
Figure 12 - for [Si/H] < -1 (Prantzos 1996, Fig. 4b; Prantzos & Silk 1998, Fig. 
3; Matteuci, Romano & Molaro 1999, Fig. 7b). The D/H declines at higher 
abundances, because a significant fraction of the atoms in the ISM have then 
come from within stars. The abundance at which the D/H decline becomes 
significant depends on the details, including the element used to gauge the metal 
abundance, the amount of that element ejected from stars, the fraction of the H 
in the ISM which has not been inside a star, the infall of gas with abundances 
nearer to the primordial values, and the amount of mixing. 

The Galactic chemical evolution models apply on average to large volumes 
of gas which are well mixed. At the epoch of observation the gas may instead be 
isolated, and the elemental abundances may have been determined by relatively 
few stars. This will allow more stochastic variation in the D/H at a specific 
metal abundance, but it will not change the general expectation that D/H will 
remain near the primordial value as long as the metal abundances are low. 
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9. New Values for the Relevant Cosmological Pcirameters 

We take the weighted mean of the D/H values from three QSOs, PKS 1937- 
1009, Q1009+2956 and HS 0105+1619, as our best estimate for the primordial 
D/H ratio. We include all three QSOs because we beheve that the dispersion in 
these D/H values is dominated by systematic errors. We weight by the quoted 
errors, which we believe indicate the relative sizes of the random errors. For the 
error on the best value of D/H we take the dispersion in the three values, divided 
by Vs. We do not use the error on the weighted mean, which is unrealistically 
small, because the dispersion on the values is larger than expected. We then 
estimate: 

• log D/H = -4.52 ± 0.06 (15% error), or 

• D/H = 3.0 ±0.4x10-^ 

This new D/H value, together with over 50 years of theoretical work and 
laboratory measurements of reaction rates, leads to the following values for 
cosmological parameters (Esposito et al 2000a,b; Buries, NoUett & Turner, 
2000): 

• 7/ = 5.6 ± 0.5 X 10-1° (from standard BBN and D/H) 

• Qbh^= 0.0205 ± 0.0018 (baryon density, in units of the critical density, 9% 
error) 

• Yp = 0.2471 ± 0.0009 (predicted mass fraction of ^He, from SBBN and 
D/H) 

• ^He/H = 1.09 ± 0.06 x 10"^ (predicted from SBBN and D/H) 

• ^Li/H = S.Stl'H X 10-1° (predicted from SBBN and D/H), 

where the error on Yp includes the error from D/H (0.0008) and from the 
calculation of Yp for a given rj (0.0004) (Lopez & Turner 1999). 

These values are generally consistent with most other measurements (Tytler 
et al. 2000; Buries, Nollett & Turner 2000; Olive, Steigman & Walker 2000). For 
example, Pagel (2000) reviews Yp measurements and concludes that "systematic 
errors up to about 0.005 are still not excluded", and that "Yp is very probably 
between 0.24 and 0.25". Of special note are ^Li and the cosmic microwave 
background (CMB) estimate of Qbh^. 

The predicted value for the ^Li/H abundance is a factor of two to five times 
higher than measured in halo stars (Bonifacio & Molaro 1997; Ryan, Norris & 
Beers 1999). Models allow at most a factor of two destruction of the ^Li in these 
stars (Ryan et al. 1999, 2000; Deliyannis & Ryan 2000; Tytler et al. 2000). 
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Recent measurements of Jl^/i^ from the CMB are in agreement with our value 
from D/H and SBBN. Compared to our more accurate Qf^h"^ value, the value 
from BOOMERANG (Lange et al. 2000: nbh'^= 0.031 ±0.004) is approximately 
3a higher, that from MAXIMA-I (Balbi et al. 2000: nbh'^= 0.025 ± 0.010) is 
0.5a higher and that from the Cosmic Background Imager (Padin et al. 2001: 
flhh^= 0.009) is just over la lower. The agreement of these flbh"^ values is a 
dramatic validation of the physics used in SBBN and CMB. 

10. Summary & Discussion 

We have presented the fourth quasar, HS 0105+1619, which shows an 
absorption system having a low deuterium to hydrogen abundance ratio: 
D/H= 2.54 ± 0.23 x 10~^ in the z ~ 2.536 Lyman limit absorption system. We 
first obtained low resolution spectra in our survey for D/H QSOs with the Lick 
3 meter Shane telescope, and here we presented over 24 hours of spectra from 
the HIRES spectrograph on the Keck-I telescope. 

The absorber has a high neutral hydrogen column density, 
log Nhi— 19.422 ± 0.009 cm~^ which is 36 times larger than the next 
highest case studied for D/H, but 8-100 times less than standard damped 
Lyman alpha systems. Very little is known about absorbers with log Nhi— 19 
cm^^, and this absorber may not be representative of this class because it was 
selected to have a very simple velocity structure and low b values. 

While the absorber towards object has by far the highest neutral H I column 
density of the absorbers which we have studied for D/H, it also has the lowest 
total Hydrogen column density, by a factor of 4 - 7, when we correct for the 
ionization. 

We know little about the environment around the absorber. It might be in 
the outer parts of galaxy, as are LLS at low rcdshift, or in a disk, as are damped 
Lyman alpha absorbers at high redshift. Alternatively, it may be in a relatively 
isolated gas cloud, because see just one component, with an exceptionally small 
spread of velocities. In either case, numerical simulations of the growth of 
structure suggest that the absorbing gas has been incorporated into a galaxy by 
today. 

For the first time, we detect deuterium absorption in 5 Lyman series 
transitions, and determine log Ndi— 14.826 ± 0.039 cm~^. We have strong 
arguments that the observed absorption is indeed deuterium, and not inter-loping 

hydrogen or metal line absorption. 

We observe a number of associated metal line absorbers, from which we 



-17- 



calculate that the gas is warm and neutral. The metaUicity of the system 
is ~ 0.01 times solar, indicating that the measured D/H is representative of 
primordial D/H. 

We argue that HS 0105+1619 offers the most secure detection of D/H to 
date, and that the D/H ratio determined from all QSOs has been made more 
secure. 

10.1. HS 0105+1619 gives the most secure mecisurement of 

primordial D/H 

The measurement of primordial D/H towards HS 0105+1619 is more 
secure than our prior measurements towards PKS 1937-1009, Q1009+2956 and 
QO 130-4021 for several reasons. By secure we mean that we have the most 
information, and hence there is less chance of undetected errors which might 
greatly exceed those quoted. We are not exphcitly referring to the size of the 
quoted errors, which are similar for PKS 1937-1009 and HS 0105+1619, and 
larger for Q1009+2956 . 

The absorption system in HS 0105+1619 is simple. Like Q0130-4021, the 
second most secure result, the absorber in HS 0105+1619 is modeled with a 
single component, which simplifies the measurement of the column densities. 
PKS 1937-1009 was modeled with two or three components, and Q1009+2956 
with two to four. 

HS 0105+1619 is the only case to show more than one strong deuterium line, 
which reduces the chance of contamination, and gives more reliable b values, 
and hence Ndi- 

For HS 0105+1619 we listed above the many reasons why the absorption near 
the deuterium line position is deuterium. The chance of serious contamination 
from the Lya forest will decrease with rising N^j. Such contamination is least 
likely to be significant in HS 0105+1619, followed by PKS 1937-1009. There 
is some contamination in Q1009+2956 while for Q0130-4021 we see a lot of 
contamination, and obtain only an upper limit on Noi. We do not know 
whether the chance of contamination by components of the Lyman limit system 
changes with Nhi- 

We also expect that the chance of contamination decreases as the b value 
of the deuterium decreases, because H lines often have larger b values. Hence 

HS 0105+1619 is the most secure detection of deuterium. 

For HS 0105+1619 we have the most information on the velocity field and b 
values because we see several D I, N I and O I lines. 
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For HS 0105+1619 the metal abundance is obtained with additional 
redundancy. The gas is nearly neutral, and hence wc get the abundances of 
several elements: H, C, N, O, Al, Fe and Si. For the PKS 1937-1009 and 
Q1009+2956 we used a standard photoionization model to find the level of 
ionization which explained the relative abundances of ions such as C II, C III, 
C IV and Si II, Si III, Si IV. We obtained a solution for each element, and 
these agreed, which provided a check. For Q0130-4021 the ionization and metal 
abundances are both less well known. 

10.2. The primordial D/H becomes more secure 

The new measurement makes the primordial D/H much more secure because 
in each case we are sampling gas with different physical conditions, and some 
systematic errors, including those associated with the measurement of column 
densities, may be different for each QSO. First, the new measurement is the 
most secure. Second, we have increased the number of QSOs in which we have 
measured deuterium from two to three. Third, we sample a new region of 
space. Each absorber is a different direction in the universe, and each samples 
a sight line of about 1-10 kpc, which requires about 10^ — 10^ solar masses of 
gas. Fourth, the absorbing gas covers a factor of 240 range in Nm. Fifth, the 
absorption systems have differing ionization, with a range of 2000 in the H I/H 
ratio. 

There is less variation in other parameters. The metal abundances cover a 
factor of ten, which is the typical for QSO absorption line systems, while the 
redshifts cover most of the range observable from the ground. 
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Fig. 1. — Spectrum of HS 0105+1619. The upper panel shows the low resolution flux 
calibrated spectrum obtained with the Kast spectrograph. The lower panel shows the flux 
calibrated HIRES spectrum. The flux calibration was noisy at wavelengths less than 3800 A 
and was not applied to the Lyman limit, which is not shown above for the HIRES spectrum. 
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Relative Velocity [km/s] 

Fig. 2. — Lyman scries absorption in the z ~ 2.536 Lyman limit system towards 
HS 0105+1619. The velocities shown are relative to the H I redshift oi z = 2.535998. 
The vertical scales are linear flux, from zero, and the lower traces are the Icr error. The D I 
absorption is seen at —82 km s~^ in Ly/3, Ly7, Ly-5, Ly-6 and Ly-7. 
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Fig. 3. — Observed metal line absorption associated with the z ~ 2.536 Lyman hmit system. 
The metals are grouped according to ionization state and are organized by atomic mass. 
The low ionization lines have simple, narrow profiles centered near km s~^. 
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Fig. 4. — Spectral regions used to measure the H I column density. The upper panel shows 
the Lya line, with the core and damping wing regions used in the fit shaded grey. The lower 
panel shows the absorption near the Lyman limit. Over-layed is the single component fit to 
the hydrogen with a column density of log Nhi= 19.422 cm~^ from Table ^ . 
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Fig. 5. — Lja region of the z ~ 2.536 Lyman limit system. Over-layed is the continuum fit 
(sohd fine) and the approximate la error to the continuum fit (dashed fine). 
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Fig. 7. — Simultaneous fit to the deuterium at —82 km in 5 Lyman series transitions 
with logA^/)/= 14.81 cm~^ and b = 9.93 km s~^. Also included in the fit is the H at km s~^ 
and additional Lya forest absorption. 
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Fig. 8. — Values of the deuterium column density, log N^i (cm~^), and the absorption width 
parameter, 6, for each of the 5 deuterium lines fit separately. The bold cross represents the 
values when all lines are fit simultaneously using the best estimate for the continuum level. 
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Fig. 9. — Best fit to metal lines associated with the z ~ 2.536 LLS towards HS 0105+1619. 
The velocities shown are relative to 2; = 2.535998. 
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Fig. 10. — Intrinsic velocity widths measured for the neutral and low ionization ions. The 
straight line is the best fit to H I, N I, and O I alone. The slope of the line gives the 
temperature of the gas, T — 1.15 x 10^ K, and the intercept gives the turbulent velocity, 
hturh = 2.56 km s~^. The D I absorption comes from the same gas as the H I, N I and O I. 
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Fig. 11. — The predicted ionization correction loge(X/H I) for ion X as a function of the 
intensity of the photoionizing radiation, log J912, in units of ergs cm~^ Hz~^ sr~^. We 
conclude that the ionization is low, as on the left of the plot. 



-33- 



(X) 



in 

I 

O 



X 

Q 



O 



T 1 1 1 1 1 1 1 1 1 1 1 1 1 1~ 

'q0130-4021 



Q1009+2956 



PKS1937-1009 



HS0105+1619 



lismI 



■3 



2 







[Si/H] 



Fig. 12. — The relationship between D/H measurements and metal abundances, as gauged 
by sihcon. Standard chemical evolution models predict httle change in D/H until the [Si/H] 
> — 1, at which time a significant fraction of the gas has come from inside stars, where the D 
has been destroyed. The values of [Si/H] for PKS 1937-1009 and Q1009+2956 represent the 
mean value, weighted by the total H column density in each component: —2.26 and —2.53 
respectively. The D/H value for the local interstellar medium (LISM) is from Linsky (1998), 
and the LISM [Si/H] is from Savage & Sembach (1996). 
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Table 1. OBSERVATIONS OF HS 0105+1619 



Instrument 


Date 


Integration Time 


Wavelengths covered 








(seconds) 


(A) 


Kast 


August 16, 


1998 


2700 


3200 


- 5100 


HIRES 


October 10 


, 1999 


7200 


3200 


- 4720 


HIRES 


October 11 


1999 


2 X 8000 


3200 


- 4720 


HIRES 


November 9, 1999 


1800 


4220 


- 6640 ^ 


HIRES 


September 


19, 2000 


2 X 7200 


3200 


- 4720 


HIRES 


September 


20, 2000 


8600, 10800 


3200 


- 4720 


HIRES 


September 


20, 2000 


4 X 7000 


3200 


- 4720 


Total (HIRES): 






86,800 







^The wavelength coverage for this observation was not continuous due to 1-10 A spectral order gaps. 
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Table 2. IONS OBSERVED IN THE z ~ 2.536 LLS TOWARDS 

HS 0105+1619^ 



Ion 


log N 




z 






(cm~^) 


(km s~^) 




(km s~^) 


H I 


19.422 ± 0.009^^ 


13.99 ±0.20 


2.535998 ± 0.000007 


0.0 ±0.6^ 


T 

D 1 


14.826 ± 0.039 


9.85 ± 0.42 


o r'o/^i^i^o 1 r\ r\r\r\r\r\Ci 

2.536002 ± 0.000008 


0.4 ± 0.7 


N I 


12.306 ±0.060 


5.13 ± 1.57 


2.535998 ± 0.000009 


0.0 ±0.8 


I 


14.378 ± 0.024 


4.30 ±0.11 


2.535991 ±0.000001 


-0.6 ±0.1 


C II 


14.349 ±0.069 


5.44 ±0.21 


2.535985 ± 0.000001 


-1.1 ±0.1 


Si II 


13.156 ±0.012 


5.01 ±0.09 


2.535980 ± 0.000001 


-1.5 ±0.1 


Al II 


11.870 ±0.040 


9.69 ±0.99 


2.536019 ± 0.000009 


±1.8 ±0.8 


Fell 


12.913 ±0.087 


6.45 ± 2.33 


2.535983 ± 0.000014 


-1.3 ± 1.2 


C III 


13.716 ± 0.032 


14.48 ±0.56 


2.535923 ± 0.000004 


-6.4 ±0.4 


C IVs 


13.277 ±0.019 


20.40 ± 1.07 


2.535913 ± 0.000009 


-7.2 ±0.8 


N II 


13.501 ±0.013 


6.82 ±0.30 


2.536015 ±0.000002 


+1.4 ±0.2 


Si III 


13.097 ±0.037 


7.10 ±0.24 


2.535977 ±0.000001 


-1.8 ±0.1 


Si IVs 


12.647 ±0.053 


9.16 ± 1.59 


2.535995 ± 0.000014 


-0.3 ± 1.2 



^ The three sections group ions by increasing ionization. For ions in the same gas, we expect 
the h values to decrease with increasing mass. Errors quoted in the table for the N and z values 
are from VPFIT alone, except for H I and D I. 

^ The intrinsic h value. 

Velocities are all relative to ^; = 2.535998. The errors in v values come from the errors 

listed on the z values: cr(w)^ = (c(t(z)/(1 + z))^ + 0.09^, where 0.09 km s"'^ is the minimum 
internal uncertainty in the wavelength scale. The internal error may be 1 - 2 km s~^, while 
the external error is approximately ±10 km s~^. 

^ The value for log Nhi is the weighted mean from the damping wings and core region of 
Lya, and the error includes the continuum level error. 

^ The HI z defines v = 0. The 0.6 km error is the uncertainty in the v of the H lines 
in this frame. 

The A'' and b for the D I are the weighted means of the individual fits to the five D 
transitions, and the error on the N includes the contribution from the continuum uncertainty. 

^ The lines of this ion may be multiple component blends. 
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Table 3. SPECTRAL REGIONS USED TO MEASURE D I 



Region 


^min (-^) 




Ly/3 


3625.1 


3627.4 


Ly7 


3437.2 


3439.2 


Ly-5 


3314.2 


3316.4 


Ly-6 


3289.4 


3291.4 


Ly-7 


3273.5 


3275.2 



Table 4. ADDITIONAL LINES USED TO FIT D REGIONS 



logA^ (cm ^) 


b (km s ^) 


z 


Region'^ 


14.324 


25.2 


2.53455 


all^ 


14.901 


18.2 


2.53554 


all'' 


13.365 


26.2 


1.72637 


Ly-5 


12.944 


29.7 


1.70682 


Ly-6 


12.770 


39.7 


1.69376 


Ly-7 



'^The regions are defined in Table || 

^All D regions used: Lya, Ly/3, Lyj, Ly-5, Ly-6, and Ly-7 
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Table 5. PARAMETERS FOR THE D/H MEASUREMENTS 



QSO 


logD/H 




log Nhi 


log uh i/uh 


m 


[Si/H] 








(cm~^) 




(km s~^) 




PKS 1937-1009 ^ 


-4.49 ± 0.04 


3.572 


17.86 ±0.02 


-2.35,-2.29 


14.0 ± 1.0 


-2.7,-1.9 


Q1009+2956 ^ 


-4.40 ±0.07 


2.504 


17.39 ±0.06 


-2.97, -2.84 


15.7 ±2.1 


-2.4,-2.7 


Q0130-4021 


< -4.17 


2.799 


16.66 ±0.02 


-3.4 


16-23 


-2.6 


HS 0105+1619 ^ 


-4.60 ±0.04 


2.536 


19.422 ± 0.009 


-0.1 


9.85 ± 0.42 


-1.85 



We list the parameters for each of the two components, where available. Results from Tytler, Fan & Buries 

(1996); Buries & Tytler (1998a). 

^ Wc list the parameters for each of the two components, where available. Results from Tytler & Buries (1997) 

and Buries & Tytler (1998b). 

^From Kirkman et al. (1999). 
^This paper. 



